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SUMMARY 

Analysis for the turbulent mixing of free jets is presented in this paper and com- 
pared to recent experimental results. A turbulent mass diffusion model is presented 
and is based on the concentration potential core. The model yielded good results when 
compared with the experimental results except for low-speed flows where few experi- 
mental data are available. 

A review of recent experimental results verifies again that the three diffusion pro- 
cesses in turbulent mixing are interrelated; however, no single diffusion model may be 
used for all three processes. This is especially true when pressure gradients are 
present in the flow field. It is shown that even though momentum diffusion is signif- 
icantly affected by pressure gradients, mass diffusion is not. 

It is further indicated that the mass diffusion model has been derived and is 
based on the accurate correlations of experimental results obtained for the concentra- 
tion potential core. Similar techniques may be used in deriving an expression for the 
momentum and thermal diffusion coefficients. These expressions would be more com- 
plicated since they would have to take care of boundary layer at the start of the mixing 
region. 

Finally, a comparison of the analyses, using this particular model and Ferri's 
model, with available experimental results is made. 

INTRODUCTION 

The problem of free turbulent mixing with large density gradients has been an 
area of considerable interest and ever increasing practical importance in the past 
decade. Problems of constant density mixing may be traced back as far as the 1930's. 
These types of problems have a wide variety of applications and have been used for 
wakes, rocket planes, supersonic combustion, nuclear core reactions, and so forth. 


* This research was sponsored in part by the Aerospace Research Laboratories, 
Air Force Systems Command, United States Air Force, Wright -Patterson Air Force 
Base, Ohio, under Contract F-33615-68-C-1184 and in part by NASA International 
Fellowships. 


In reference 1, the Prandtl eddy -viscosity model for constant density flow was 
extended to the variable density flow. However, as pointed out in reference 2, this model 
fails when the mass flux of each stream is equal. Ting and Libby (ref. 3) determined the 
compressibility factor for the Prandtl model using their transformation, while Donaldson 
and Gray (ref. 4) and Peters (ref. 5) obtained the compressibility factor empirically. 

These models, however, fail also when the velocity of each stream is equal, since the 
Prandtl model becomes zero. Zakkay, Krause, and Woo (ref. 6) developed the eddy- 
viscosity model in terms of the center-line velocity and velocity half-radius. Schetz 
(ref. 7) extended the Clauser model (ref. 8) to the axisymmetric wake which was 
expressed in terms of mass defect across the wake. 

Each of all these models has a constant which is chosen such that the numerical 
solutions of wake equations agree with the experimental data by employing the proposed 
eddy viscosity and appropriate constant values for turbulent Prandtl number Pr^ and 
turbulent Schmidt number Sc t> Therefore, it has been found in references 2,9, and 10 
and also in the present investigation that all these models provide good predictions for 
certain very restricted flow conditions, but none of them is valid for the general case. 

In comparison with these numerous works for the eddy viscosity, only few exten- 
sive works have been done for other transport properties, namely Pr t , Sc t , and Le t 
(turbulent Lewis number). In analytical and numerical investigations so far, constant 
values for Pr t and Sc t have been assumed, and several experimental studies seem 
to approve these assumptions of constant properties under the condition of zero pressure 
gradient in the mixing flow field. 

In reference 11, Forstall and Shapiro determined Le t = 1.0 and Pr^. = Sc^. * 0.7 
using the data of air-to-air mixing measured at very low speed. Zakkay et al. (ref. 6) 
carried out extensive measurements for the turbulent mixing of high-speed coaxial jets 
comprised of several dissimilar gases. For hydrogen and air mixing, Sc t = 0.8 to 2.0 
and Let = 0.9 to 1.2 have been concluded. Chriss (ref. 12) measured very detailed pro- 
files in mixing flows of high-speed coaxial air and hydrogen jets and obtained Let “ 

The accuracy of these experimental data has been checked by Zelazny, Morgenthaler, and 
Herendeen (ref. 9) and Harsha (ref. 10) by means of constant-momentum integral check 
and was found to be very good. Peters, Chriss, and Paulk (ref. 13) have shown for one 
case of hydrogen and air mixing measured by Chriss that Prt = Set = 0.85 is approxi- 
mately valid. Zelazny et al. (ref. 9) also calculated Set f° r f° ur cases of Chriss’ data 
(ref. 12) and confirmed the result obtained by Peters et al. (ref. 13). 

It is difficult at this time to generalize and to derive equations for a single applica- 
tion; however, some basic conclusions have been reached from past research efforts. 

For example, the authors have indicated previously that the mixing which occurs between 
two coaxial streams depends largely on their initial energy and their transfer phenomena 
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in the mixing region. Therefore, differentiation must be made between jets where momen- 
tum diffusion is the controlling mechanism and jets where mass diffusion, or thermal 
diffusion, is the controlling mechanism. Theoretical treatment of the problem is impos- 
sible since turbulent transport properties for the three diffusion processes are not yet 
available. Reference 6 has attempted to deduce such turbulent transport coefficients 
from mass diffusion experiments. 

Similar attempts have been made in wakes to deduce momentum diffusion coeffi- 
cients; however, no single model has been capable of analyzing the entire range of flow 
fields. In addition to the complications of the unknown transport properties, in several 
problems, the pressure field is not constant, and therefore complicates the analysis 
further. For instance, references 14 and 15 have indicated that pressure gradients do 
not affect mass diffusion; however, momentum diffusion is significantly affected. These 
regions of varying pressure gradients are present close to the exit of the jets and are 
directly responsible for the discrepancies among constant -pressure analyses. Results 
of mixing with pressure gradients, in reference 15, indicate a large discrepancy between 
theory and experiments for momentum diffusion, although good agreement exists for mass 
diffusion. 

Therefore, since the topic under discussion is complex and cannot be described in 
terms of one single model, the present paper will concentrate on the mixing of coaxial 
jets with special emphasis on mass -diffusion processes. At first a description of the 
model used for mass diffusion will be derived (see fig. 1) and then the predictions of this 
model are compared with several of the test cases supplied by the organizers of the 
Langley working conference on free turbulent shear flows. 

SYMBOLS 

Cp specific heat at constant pressure 

Cp normalized specific heat with respect to c p ^ 

®Pj^ normalized specific heat of kth gas; k = 1 corresponds to inner jet gas 

and k = 2 corresponds to outer jet gas 

Dj inner jet nozzle diameter 

Dt turbulent diffusion coefficient 

H total enthalpy of gas 
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M 
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Prt 

R 

Re 

r 

r 

r i 

r mc 

%c 

r mu 

Sc t 

T 

T 

u 
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static enthalpy of kth gas 
turbulent Lewis number 
local Mach number 

total number of grid points in T-direction 
concentration decay exponent along center line 
static pressure 
turbulent Prandtl number 
gas constant 
Reynolds number 
radial coordinate 

normalized radial coordinate with respect to 

radius of inner jet nozzle 

concentration half -radius defined so that Yj 

normalized half -radius with respect to rj 

u* + u e 

velocity half-radius defined so that u = - at r = r mu 

turbulent Schmidt number 
temperature 

normalized temperature with respect to T e 

axial component of velocity 

normalized axial velocity with respect to u e 


* 

= 2 at r = r mc 
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radial component of velocity 


nondimens ional radial velocity 
ratio of molecular weights, Wjy/W e 
normalized molecular weight with respect to W e 
axial coordinate 

initial value of x to start computation 
concentration potential core length 
normalized axial coordinate with respect to rj 
normalized x Q with respect to r^ 
normalized axial coordinate with respect to x 0 
mass concentration of inner jet gas 

mass concentration of kth gas; k = 1 corresponds to inner jet gas, and 
k = 2 corresponds to outer jet gas 


transformed axial coordinate defined as 



die 


eddy viscosity 


specific -heat ratio 


mass -flow ratio defined as 


P i u i 


Pe u e 


modified stream function defined by equation (7) 


normalized center -line velocity defined by 
transformed axial coordinate defined as 2 


u e - u<t 




density 


H 


(pDt) 


normalized density with respect to p e 
modified stream function defined as 
turbulent shear stress 


is: 


\ 0.5 


pur dr 


turbulent viscosity 

normalized value of pD t with respect to pjUjrj 

normalized value of pu with respect to p.Uj 

J J 


(pu) 

(p V ) normalized value of pv with respect to p g u e 

Subscripts: 


t center-line values 

e external conditions 

i initial conditions 

j inner jet conditions 


THEORETICAL ANALYSIS 


Basic Equations and Turbulent Diffusion Coefficient 

The governing equations for the mean turbulent flow properties in coaxial mixing 
with axial pressure gradient can be described as follows: 

Conservation of mass: 

!^ PU) + r lr' pVV) = 0 (1) 

Conservation of momentum: 


9u 9u 
PU ^ + PV ^ = 



dp 

dx 


( 2 ) 
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Conservation of energy: 


9H 9H , dp 13 pe _ 3H 
' ,u '5r + pv 9F = u te + ?-5?|p^ r 'to 




Pr. - 1 


+ ^T peru ^ + 


9u 

9r 


Le t - 1 
Pr t 


per 




(3) 


Conservation of species: 


9Yi 


pu _i +pv !Ii = i(£ir^i] 

P 3x P 9r r \Sc t 9r ) 


( 4 ) 


where the following conditions were assumed: 


(1) The flow is chemically frozen. 

(2) The radial pressure gradient is not existing. 

(3) The streamwise pressure gradient is supposed to exist and therefore the 
external flow conditions could be variable. 

In references 16 and 17, equations (1) to (4) are solved for the axisymmetric coaxial 
jet mixing with the assumptions of dp/dx - 0, Pr t = Sc t = Le t = 1.0, and step initial pro- 
files for velocity, enthalpy, and concentration. Obviously these assumptions reduce those 
equations to a single equation and the solution can be applicable for any flow properties. 
As for these assumptions, however, the transport coefficients Pr t , Sc t , and Le t are 
generally not equal to unity, and the initial profiles are not always step profiles because 
of boundary layers developed on the jet nozzle wall. The exact application of this solu- 
tion is possible only for the species equation involving the turbulent diffusion coefficient, 
for which the initial profile can be assumed correctly to be a step profile and in which 
other transport properties are not directly involved. Furthermore, even if the axial 
pressure gradients exist, the mass diffusion will not be affected significantly as shown 
experimentally in references 6, 14, and 15 and in the present experiment, since the pres- 
sure gradient affects the concentration indirectly through velocity and temperature field 
as seen in equation (4). Therefore, the solutions as obtained in references 16 and 17 can 
be employed exactly only for the species equation involving the turbulent diffusion coeffi- 
cient in the mixing with axial pressure gradients. 

The species equation with the diffusion coefficient is 


9Yi 1 9 / 


( 5 ) 


where the initial conditions at x = 0 are 
v. = 1 


(for 0 = r a r<'i 


Y i=° 

and the boundary conditions are 

9Y i 

— ±= 0 

dr 

O 

A 

A 

(at r = 0) 

o 

II 

(at r — <») 


Applying the modified von Mises transformation, Kleinstein obtained the exact solu- 
tion for the linearized boundary -layer equation in reference 17. Employing the same 
method, the species equation can be solved as follows: 


Equation (5) is transformed to 


d J± m ±J_ 

dx 'I' 


( P D,)(pu)f 2 (l^ 


by the transformation terms: 







( 6 ) 


(7) 


Neglecting the fourth -order term in the Taylor series expansion of 


(pDt) Y * °(* 4 ) 

the species equation (6) reduces to the linear form in -plane as 


* 8 *\ wj 


where £ is defined as 



Assuming a step initial profile, the solution of this well-known heat -conduction 
equation can be expressed along the center line by 


Yj =l-exp - 



( 8 ) 


where \kj is \[2 by the definition given in equation (7). In the far wake region as 4 — 
approaches the Gaussian distribution (ref. 6) as 
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= exp \'4T, 


These solutions (eqs. (8) and (9)) can be used to introduce the diffusion model in a 
way similar to that developed in reference 6 where Zakkay et al. introduced the eddy- 
viscosity model. 

Since (pu) — i in the far wake region as t, — 

* 2 _ 2f * 

* ~T~ 


= exp 


Therefore, assuming the center-line concentration decay law to be Yj^ = x n = , 

the concentration profile can be expressed in the physical plane by 

Yj = *- n (l - x- n ) r /X (11) 

Then the asymptotic jet spread can be given by taking the limit of equation (11) as follows: 

Um *4x ln lLY /2 | x > 1 ) 

Integrating this equation gives the half -radius for concentration r mc in the far wake 
region as 


-(*) 


r mc - 


nll/2 


Equation (12) may provide an analytical prediction of the concentration half-radius. Once 
the concentration profile is expressed in the physical plane by equation (11), the turbulent 
diffusion coefficient can be derived from the species equation inversely. 

The species equation (5) is described with the aid of the continuity equation by 

Hi 4 r M r .dr' -Y, f r'^ldr' 

* dr J *^0 dK 

Then the diffusion coefficient along the center line can be expressed as 




'a. S 


9 Yi 
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Substituting equation (11) into equation (13) gives the turbulent diffusion coefficient on the 
center line as follows: 

n-1 

<u, 

In equation (14) there are two unknown parameters, namely x Q and n. The decay 
exponent n in Yj^ = x has been shbwn to be 

n = 2 (15) 

for high-speed jet mixing in many experiments (for example, refs. 6, 14, and 15) and also 
will be shown to be n = 2 in the present experiment. 

Now the potential core length Xq for concentration in high-speed mixing, where 
n = 2 is valid, should be formulated in terms of known values. The parametric inves- 
tigation has shown the importance of momentum ratio at the initial exit plane on the poten- 
tial core length for concentration. The typical correlation between the potential core 
length for concentration and the initial Mach number ratio can be seen obviously in fig- 
ure 2. The Mach number ratio of two jets can be expressed by the momentum ratio as 

PjUj 2 _ (if)] /Mj\2 

P e u e 2 fafi ~W 

\RT/e 

since Pj = p e . Therefore, figure 2 shows that when the momentum ratio becomes larger, 
the potential core length becomes longer. 

These data are rearranged in figure 3, from which the correlation equation has been 
obtained as 


ef) 

/concentration 


= 1 + 301-4- 


1/4 


Mr 


1 + Mj exp | 



(16) 


From equation (16), when Mi = 0, 

\ ^/concentration 

which corresponds to the axisymmetric wake, and the result of x 0/ /t>j = 1.0 is supported 
by the data of the supersonic wake studies (refs. 18, 19, and 20). When M e = 0, which 
corresponds to the jet injected into the quiescent atmosphere, the data of Keagy and 
Weller (ref. 21) and O'Connor (ref. 22) are fairly well correlated as seen in this figure. 

Now the turbulent diffusion coefficient, assuming that it is only a function of x, can 
be described with the aid of equations (14) and (15) as 
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It is noticeable that equation (17) does not include any adjustable constant to match the 
numerical solutions of wake equations with the experimental data. Such constants were 
included in all previous viscosity models suggested by Prandtl, Ferri, Schetz, and others. 


As shown subsequently the diffusion model given by equation (17) provides good 
numerical predictions for the high-speed mixing of hydrogen and air, for which the decay 
exponent n satisfied equation (15) approximately. However, for the general case of 
arbitrary dissimilar gases, it has been found necessary (by means of numerical calcula- 
tions) to include a molecular -weight factor given as 1 + 0.8(Wj/W e ) in the expression 
for turbulent diffusivity; thus, 

p Dt = °- 6 - £!i 

1 1 + 0.8W jL 


me 


(Pu) d 


(18) 


.1 


where 


*o = 2 


/Pi V 
+ 60-U 

\r a / 


1/4 


Mf 2 1 1 + Mf exp 


J 


3 


Me \ 2 

Mi 


n~V2 




Since the ratio of the molecular weights W = Wj/W e is very small for the mixing of 
light gas and heavy gas, equation (18) becomes of the same form as equation (17) as, for 
example, in the case of hydrogen-air mixing. 


It is important to realize that the model chosen is based on the potential core 
derived from a step profile. Therefore, one cannot expect it to work if the initial pro- 
file consists of a boundary -layer profile. However, these difficulties could be avoided 
by always starting with an initial step profile, and calculating the length of the potential 
core for the step profile. This distance should be added to the calculations performed 
with the boundary -layer type profile. This procedure will result in a shift of the curve, 
which will indicate a slower mixing region. This fact will be evident in some of the 
results that will be given in the next section. 


Application to the Numerical Solution of Wake Equations 
Without Pressure Gradients 

In order to investigate the validity of the present diffusion model, the wake equa- 
tions have been solved for the cases without pressure gradients. The wake equations 

1 For numerical calculation, the concentration half -radius r mc should be deter- 
mined by the numerical procedure instead of equation (12), which will provide the ana- 
lytical prediction of the half-radius in the far wake region only. 



in the general form are written in the modified von Mises plane as follows: 


Momentum equation: 

/- - 2 ' 
9u _ 1 9 PP t Ur 9u 

\ u 9u e 

9£ 7 97 \ 7 97/ 

/"u e 9? 

Energy equation: 

/ 


9T_ 1 8%^ 

! 4 +( > 

9 ? 7Cp 9^^ Pr t 7 

dT J V 

2_ 

| p ^t r u 9T V 

c 9Yk 

SCfC r2 9 t Lj 

to k 

P k 9 t 
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T - 1 9u e PM t r U /9u 


pcp u e 9 ? c t 2 V 9t . 


Species equation: 

9Y k_ 1 3 ( p ^ 2 9Y k) 
9£ T 9 t\ Sc.t Bt . 


(k =1,2,. . .) 


where the modified von Mises transformation is defined by 


97^ 

— — = pur 
9r 


3r^ 

3x 


= -pvr 




Re~* dx 


and = Sc^pD^ and Pr^. = Le^Scj. are used. 

The associated initial and boundary conditions are 
u(0,t) = u^t) 

T(0,t) = Ti(r) 

Yj(0,7) = Yji( T ) 


^u«,0).^T(?,0).^Y,(?,0).0 
lim u(^,r) = lim T(C,r) = 1.0 

7-*-°0 7"_»oO 


9 T7T 


(19) 


( 20 ) 


( 21 ) 


lim Yj(^,t) = 0 

7-»oo J 
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The numerical calculations for equations (19) to (21) were carried out by an implicit 
finite -difference scheme. The details of this numerical process have been reported in 
references 23 and 24. Also, some pertinent results relating to this program have been 
reported in references 25 to 27. 

Once initial and boundary conditions are specified, equations (19) to (21) can be 
solved with the aid of equation (18) for diffusion coefficient and with the assumptions of 
appropriate constant values for Sc t and Pr t . In this present calculation, 

Sc t = Pr t = 0.8 and Le t = 1.0 have been used. 

The numerical calculations were carried out for the experimental data of Chriss 
(ref. 12) and Alpinieri (ref. 2). In figure 4, the numerical solutions for the data of Chriss 
are shown compared with the experimental data which were measured in hydrogen and air 
mixing. The accuracy of those experiments have been found to be very good by Zelazny 
et al. (ref. 9) and Harsha (ref. 10) (reportedly within 4 percent). Since the lateral profiles 
measured at several stations of x were reported in reference 12, the initial profiles for 
this numerical calculation were chosen at the position of x A indicated in figures instead 
of the assumed profiles at the jet exit plane. 

Figures 4(a), 4(c), and 4(f) show very good agreement of numerical solutions with 
the experiment for concentration decay on the center line and for the concentration half- 
radius. Figure 4(e) shows an example for the lateral concentration profiles which agree 
quite well with the experimental data. Also in figures 4(a), 4(c), and 4(f) the analytical 
predictions of concentration half -radius r mc given by equation (12) are shown. The 
agreement with the experimental data is excellent, considering that equation (12) does not 
include any adjustable constant to match the solutions with data. 

Concerning the velocity profiles for the data of Chriss (ref. 12), the center -line pro- 
file and the half -radius have been predicted well as seen in figures 4(b), 4(d), and 4(g). 
These results prove that if there is no pressure gradient in the flow field the assumptions 
of Sc t = Pr t = 0.8 and Let = 1.0 are valid for the high-speed mixing of hydrogen and 
air. 

In figures 4(c) and 4(f), the numerical solutions obtained by employing the eddy- 
viscosity model of reference 1 are presented. The comparisons with the data show that 
the Ferri model provides fairly good predictions in the far wake region. 

The numerical result for the data of Alpinieri (ref. 2), which are characterized by 
the very small momentum of the inner hydrogen jet, is shown in figure 5. The result 
given by the present diffusion model agrees fairly well with the data. 

In figure 6 , the results are shown for the case of C0 2 and air mixing measured by 
Alpinieri. In these cases numerical results agree fairly well with the data which have the 
decay exponent n ~ 1.5. 
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Throughout these numerical computations the momentum integral of the wake has 
been checked at each station of finite -difference calculation and found to be constant with 
a variation of less than ±1 percent. 

COMPARISONS WITH SELECTED CASES 

In addition to the cases presented in the previous section of this report, six cases 
(identified as 9, 10, 11, 12, 20, and 21) will be chosen for this section. Each of these cases 
will be discussed separately. An additional case was chosen in order to evaluate the 
numerical techniques used in the analysis. Table I presents initial conditions for all the 
above-mentioned cases. 


Discussion of the Results of Selected Cases 

Additional case .- For the additional case (see page 735 for explanation), a step pro- 
file has been given for the velocity as well as the concentration. The velocity did not pre - 
sent any problem; however, the concentration profile had to be rounded off slightly in order 
to start the calculation. For this case, the transformed r coordinate, which is indicated 
as t, had a step size of At equal to 0.009. The total number of grid points in the 
r-direction was 60. The total time required for the calculation on the 6600 computer at 
NYU was 180 seconds. The results of this case are shown in figure 7. 

Case 9 .- In case 9, the analysis was carried out for the experimental conditions of 
Forstall and Shapiro (ref. 11). Initial profiles were given both for the velocity and con- 
centration. However, one could hardly consider this case as the mixing of two dissimilar 
gases since the center jet contained 1^- percent, by mass, of a helium tracer gas. 

Ci 

The momentum controls the mixing in this case, and therefore it dominates the mix- 
ing primarily. The results of the present analysis are presented in figure 8 and compared 
with the experimental results. It is noticed that the Ferri model, which essentially reduces 
to the Prandtl model for low speeds, has the best agreement with the experimental results. 

Case 10 .- In case 10 (see ref. 12) the initial profiles were given, and therefore no 
adjustments were needed. Figure 9 presents results of the Ferri model as well as the 
results of the present model. Both models seem to agree initially; however, the present 
model seems to give the correct trend of the experimental results. For this case, even 
though the flow is subsonic, the analysis still agrees excellently with the experimental 
results. 

Case 11 .- The results for case 11 are presented in figure 10. The center jet in this 
case is a subsonic air stream traveling at Mach 0.90 surrounded by an annular Mach 1.3 
nozzle. In this case, a 1 -percent ethylene tracer gas was mixed into the central stream. 
As far as the mixing is concerned, this was a predominantly momentum transfer, and 
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tlerefore the analysis which is presented herein and which has been derived from a con- 
centration potential core would hardly be suitable. However, the analysis was still carried 
out for this model. Two different initial conditions were chosen for case 11: 

(a) Assuming an initial step profile 

(b) Assuming the initial profile provided 

It is clearly seen that condition (a) provides the fastest mixing, and the boundary - 
layer type profile (condition (b)) results in a slower mixing region. In order to use this 
model for the potential core, the length required to reach the boundary -layer profile should 
be added as initial conditions, since the analysis here assumes a step profile. The results 
of such a correction are included in figure 10 and labeled "New model with correction." 
Even though the analysis of such a technique improves the agreement, it still falls short 
of the experimental results. 

Case 12.- The results for case 12 are presented in figure 11. The center jet is 
entirely H 2 at a Mach number of 0.89. As seen in figure 11, the present analysis agrees 
fairly well with the experimental results; however, the Ferri model falls quite short of 
the experimental results. 

Case 20.- In case 20, only a 2 percent by volume of H 2 was placed in the center jet, 
and therefore the mixing is dominated by momentum. The results for this case are pre- 
sented in figure 12, which shows no agreement with the experimental results. 

Case 21.- The center jet is entirely H 2 for case 21, and both inside and outside 
streams are subsonic. The results for this case are presented in figure 13, and the anal- 
ysis, as expected, is in fair agreement with the experimental results. 

Summary of Results of Selected Cases 

The following general results were noted for the cases selected: 

(1) For all the cases where a lighter molecular gas such as H 2 or He is placed in 
the center, the present analysis works very well (that is, X < 1). This is independent of 
the conditions whether the flow is subsonic or supersonic. 

(2) For low-speed flows, where X > 1 (predominantly flows of air-to-air mixing), 
the Ferri model seems to work best. 

(3) As noted previously, the decay such as exists for the flow described in result (1) 
is always proportional to x - ^. 

CONCLUSIONS 

A model has been presented for the prediction of the concentration decay for high- 
speed coaxial mixing. The model has been derived from accurate correlation of data 
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derived from the length of the concentration potential core. Comparisons of this model 
with available experimental results seem to indicate good agreement, except for cases 
where the mass -flow ratio is greater than 1 (mixing of gases having approximately the 
same molecular weight). For all the cases where the gas in the center was much lighter 
than the coflowing outer gas, the analyses predicted the concentration decays excellently. 
It is indicated in this paper that no single model could solve the problem of mixing, and 
attempts should be made to analyze each problem separately. In addition, differentiation 
should be made between momentum, mass, and thermal diffusion. These effects become 
much more important when pressure gradients are present in the flow field. 
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TABLE I.- INITIAL CONDITIONS FOR SELECTED CASES 


Test case 

Jet 

x i 

u j/u e 

At 

(*) 

N 

Profile 

X 

M e 

Additional 

h 2 

0.0 

3.84 


m 

u-step 

at -step (modified) 

0.262 


9 

He and air 
mixture 

0.0 




As given 

3.68 

0.0266 


H 2 

5.932 

6.3 

0.03 

60 

As given 

0.508 

0.422 

11 

Air 

0.0 

0.735 

0.03 

55 

As given 

0.634 


12 

H 2 

0.0 

2.726 

0.03 

70 

As given 



20 

H 2 and air 
mixture 

0.4568 

2.1926 


55 

As given 

2.06 


21 

«2 

5.15 

3.078 



As given 

0.398 

0.497 


* 

At step size in T-direction; t transformed f coordinate as defined by 
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Figure 1.- Turbulent jet mixing flow model. 











Figure 4.- Comparison of numerical results with the data of Chriss (ref. 12) 







Figure 4.- Continued. 
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Figure 4.- Continued. 






Figure 4.- Concluded. 
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Figure 5.- Comparison of numerical results with the Figure 6.- Comparison of numerical results with the 

data of Alpinieri (ref. 2) for the mixing of hydrogen data of Alpinieri (ref. 2) for the mixing of carbon 

and air. dioxide and air. 
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Figure 9.- Case 10. 




Figure 10.- Case 11. 








Figure 11.- Case 12. 







